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A class of calculable global models for neutrino oscillations based on Lorentz and CPT violation 
is presented. One simple example matches established neutrino data from accelerator, atmospheric, 
reactor, and solar experiments, using only two degrees of freedom instead of the usual five. A third 
degree of freedom appears in the model, and it naturally generates the MiniBooNE low-energy 
anomalies. More involved models in this class can also accommodate the LSND anomaly and 
neutrino-antineutrino differences of the MINOS type. The models predict some striking signals in 
various ongoing and future experiments. 
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I. INTRODUCTION 

The minimal Standard Model (SM) of particle physics 
contains three flavors of massless left-handed neutrinos. 
However, experiments with solar, reactor, accelerator, 
and atmospheric neutrinos have convincingly demon- 
strated the existence of neutrino flavor oscillations. This 
effect cannot be accommodated within the SM and so 
represents forceful evidence for new physics. 

A popular hypothesis attributes neutrino oscillations 
to the existence of a tiny neutrino mass matrix with off- 
diagonal components. Extending the SM to incorporate 
this notion produces a model with three flavors of massive 
neutrinos (3fSM), in which oscillations are controlled by 
a 3x 3 matrix involving six parameters: two mass-squared 
differences Atoq, Am^ tm , three angles $23, #13, and 
a phase S controlling CP violation. The first four of these 
parameters must be nonzero to match established experi- 
mental data, while recent results provide indications that 
the angle #13 must also be nonzero [1, 2]. 

In this work, we explore an alternative hypothesis 
attributing part of the observed neutrino oscillations 
to tiny Lorentz and CPT violation, which might arise 
in a Planck-scale theory unifying gravity and quantum 
physics such as string theory [3]. One motivation for 
studying alternative hypotheses for neutrino oscillations 
is based on existing data. Several neutrino experiments 
have reported potential evidence for anomalous neutrino 
oscillations that is incompatible with the 3z^SM. This in- 
cludes the LSND signal [4], the MiniBooNE low-energy 
excess [5], and neutrino-antineutrino differences in the 
MiniBooNE [6] and MINOS [7] experiments. Another 
motivation is philosophical: having more than one viable 
hypothesis is known to be of great value in guiding ex- 
perimental and theoretical investigations of new physics. 
Lorentz and CPT violation is interesting in this context 
because it naturally generates neutrino oscillations and 
moreover leads to simple global models describing all es- 
tablished and anomalous neutrino data [8, 9]. 

An appropriate theoretical framework for studying re- 
alistic signals of Lorentz violation is effective field the- 
ory [10]. In this context, CPT violation is necessarily 
accompanied by Lorentz violation [11], and the com- 
prehensive description for Lorentz and CPT violation 



containing the SM and General Relativity is given by 
the Standard- Model Extension (SME) [12, 13]. In the 
SME action, each Lorentz- violating term is a coordinate- 
independent quantity constructed from the product of a 
Lorentz- violating operator and a controlling coefficient. 
The combination of observer coordinate invariance and 
Lorentz violation implies particles in the SME follow tra- 
jectories in a pseudo-Riemann-Finsler geometry [14]. 

Over the last decade or so, many experimental analyses 
using a broad variety of techniques have been performed 
to seek nonzero SME coefficients for Lorentz and CPT 
violation [15]. The intcrfcromctric nature of particle 
oscillations suggests that sensitive neutrino or neutral- 
meson experiments might well yield the first detectable 
signals of tiny Lorentz violation. In the neutrino sector, 
recent SME-based phenomenological studies [8, 9, 16- 
33] and methodologies for experimental analysis [34, 35] 
have spurred searches for Lorentz and CPT violation by 
the LSND [36], Super-Kamiokande (SK) [37], MINOS 
[38, 39], MiniBooNE [40], and IceCube collaborations 
[41]. Searches have also been performed with neutral 
mesons [42, 43], and recent DO results suggest some evi- 
dence for anomalous CP violation [44] that could be at- 
tributed to Lorentz and CPT violation [45]. 

Here, we focus on a special class of 'puma' models in 
which the 3x3 effective hamiltonian h^ s governing oscil- 
lations of three flavors of active left-handed neutrinos is 
characterized by two simple properties: isotropic Lorentz 
violation, and a zero eigenvalue [9]. The isotropic Lorentz 
violation implies boost invariance is broken while leaving 
rotations unaffected, so h^ s is independent of the direc- 
tion of the neutrino momentum but must contain uncon- 
ventional dependence on the neutrino energy E. This 
leads to unconventional energy dependences even in vac- 
uum oscillations, producing a broad range of unique neu- 
trino behavior. The zero eigenvalue can be attributed to 
a discrete symmetry of h^ s . It ensures quadratic calcula- 
bility of the mixing matrix and of oscillation probabilities 
for all models, even when matter effects are included. 
These two features differ qualitatively from the 3^SM, 
in which the Lorentz- invariant mass terms force a 1/E 
energy dependence of all terms in h^g and the lack of 
symmetry results in calculational complexity. 

The unconventional energy dependence in h v cS gener- 
ically takes the form of polynomials in E arising from 
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Lorentz- violating operators of arbitrary dimension in the 
SME Lagrange density [46]. The polynomial coefficients 
are therefore determined in terms of SME coefficients for 
Lorentz violation. For much of this work we make the 
plausible assumption that a few terms of comparatively 
low mass dimension dominate the neutrino behavior, ci- 
ther by chance or due to the presently unknown structure 
of the underlying theory, and hence that only a few co- 
efficients are needed to reproduce the bulk of existing 
neutrino data. Indeed, the basic puma models consid- 
ered below have only three degrees of freedom, which 
includes one mass and two Lorentz-violating coefficients. 
Remarkably, two of these degrees of freedom suffice to 
reproduce all established neutrino behavior, a frugal re- 
sult compared to the five degrees of freedom required by 
the 3^SM. Moreover, the third degree of freedom nat- 
urally reproduces the anomalous results found by Mini- 
BooNE [5, 6] without introducing new particles or forces. 
Comparatively minor modifications of these simple puma 
models that preserve the discrete symmetry of h v eS can 
also accommodate the LSND signal [4] and anomalies of 
the MINOS type [7]. 

The structure of this paper is as follows. The basic 
properties of the general puma models are presented in 
Sec. II. Applications to existing experiments are dis- 
cussed in Sec. III. A specific model involving one mass 
parameter and two Lorentz-violating operators, one of 
which is CPT odd, is used for illustrative purposes. Pre- 
dictions for future experiments are presented in Sec. IV. 
Some of these are strikingly different from models based 
on the 3fSM. Variant puma models using three differ- 
ent degrees of freedom or more than three parameters 
are considered in Sec. V. Finally, Sec. VI contains some 
comments on the general nature of the models. 

The notation adopted here is that of Refs. [8, 9]. A 
mass parameter is denoted m, a coefficient for isotropic 
CPT-odd Lorentz violation is denoted , and a coeffi- 
cient for isotropic CPT-even Lorentz violation is denoted 
c( d ) , where d is the dimension of the corresponding opera- 
tor. To identify the various specific puma models accord- 
ing to their coefficient content, we introduce a convenient 
nomenclature listing coefficients in descending order of 
operator mass dimension. For example, a model with 
three degrees of freedom including a mass term m and 
coefficients and for Lorentz violation is called a 
c%a$m model. 



II. GENERAL MODEL 

In the general puma model, the effective 3x3 hamil- 
tonian h^ s describing the oscillation of three active neu- 
trino flavors e, /i, r takes the form [9] 
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where A(E), B(E), and C(E) are real functions of the 
neutrino energy E. In this work, the function A is cho- 
sen to be A = m 2 /2E 1 where m is the unique neutrino 
mass parameter in the theory. The functions B and C 
have nonstandard energy dependence, which here is taken 
to arise from Lorentz-violating terms in the SME, some 
of which may lie in the nonrenormalizable sector. The 
treatment of possible contributions to from Lorentz- 
invariant operators lies outside our present scope and will 
be given elsewhere. We assume all SME coefficients con- 
tributing to h v cS are spacetime constants, so the model (I) 
incorporates translation invariance and conserves energy 
and momentum. In the context of spontaneous Lorentz 
violation, where the SME coefficients can be interpreted 
in terms of expectation values in an underlying theory, 
this assumption implies soliton solutions, massive modes, 
and Nambu-Goldstone modes [47] are disregarded. The 
latter may play the role of the graviton [48] , the photon in 
Einstein-Maxwell theory [49] , or various new forces [50] . 
For simplicity in most specific models considered here, B 
and C are taken to be monomials in E, although more 
complicated polynomials or nonpolynomial functions can 
also be of interest. 



The function A decreases inversely with energy, while 
B and C typically increase. At low energies, the effective 
hamiltonian h v cS is therefore well approximated by the A 
term alone. This term has a 'democratic' form, exhibit- 
ing symmetry under the permutation group S3 acting 
on the three neutrino flavors e, /i, r. In contrast, the 
nonstandard energy dependences in the B and C terms 
dominate at high energies. The flavor-space structure of 
these terms breaks the S3 symmetry to its S2 subgroup 
in the \i-t sector. 

For antineutrinos, oscillations are governed by the 
CPT image of the effective hamiltonian h^g. The 
effect of the CPT transformation on is to change the 
signs of any coefficients for Lorentz violation that are 
associated with CPT-odd operators in the SME. Since 
mass terms are invariant under CPT [II], the A term in 
/ig ff is unaffected by the transformation. At low ener- 
gies, the full permutation symmetry of the puma model 
is therefore S3 x S3, where S3 is the symmetry acting on 
antineutrino flavors. At high energies, the S3 x S3 invari- 
ance breaks to S2 x S2. If any coefficients for CPT-odd 
Lorentz violation are present, differences between neutri- 
nos and antineutrinos can become manifest. 



An elegant feature of the puma model is the existence 
of a zero eigenvalue for the effective hamiltonian, which 
is a consequence of the permutation symmetry of the tex- 
ture (I). This implies considerable calculational simplifi- 
cation compared to the 3^SM and typical other neutrino- 
oscillation models. Many results can be obtained exactly 
by hand even when all three neutrino flavors mix. A short 
calculation reveals that the eigenvalues A a /, a' = 1,2,3, 
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of the effective hamiltonian h^ s take the exact form 
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(2) 



The mixing matrix U a i a that diagonalizes h" s can also 
be expressed exactly as 
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1 1 this equation, the index a ranges over a = 
the normalization factors are 

Ni = ^(\ 1 -2A) 2 + 2(A + B) 2 1 
iV 2 = v/(A 2 - 2A) 2 + 2{A + B) 2 . 



(3) 



and 



(4) 



The eigenvalues A a <, the mixing matrix U a 'a, and the 
normalization factors N\, iV 2 for the antineutrino effec- 
tive hamiltonian h v cS are obtained by CPT conjugation 
of B and C. 

In the low-energy limit, the mixing matrix (3) reduces 
to the tribimaximal form originally postulated on phc- 
nomenological grounds by Harrison, Perkins, and Scott 
[51]. The democratic structure of the A term in h^ s 
therefore ensures tribimaximal mixing of the three neu- 
trino flavors at low energies. Combined with the choice 
A = m 2 /2E > 0, this mixing guarantees agreement of 
the puma model with low-energy solar neutrinos [52] and 
with the mixing observed in KamLAND [53] . For a suit- 
able choice of mass parameter m, as discussed in the next 
section, the A term can also correctly describe the L/E 
oscillation signature observed by KamLAND [54]. 

Another defining feature of the puma model is a 
Lorentz-violating seesaw [8] that mimics a mass term at 
high energies, without invoking mass. This differs from 
the usual seesaw mechanism [55, 56], which is based on 
mass terms in the action. Suppose B and C are mono- 
mials of the form 



B{E) = k^E p - 3 , C{E) = c {q) E q - 3 , 



(5) 



where p and q are the dimensions of the operators asso- 
ciated with the coefficients k^ and & q \ In this work, 
we take > for definiteness but consider both sign 
options for k^ p \ Reversing the sign of produces phe- 
nomenology closely related to reversing instead the sign 
of k( p \ as can be seen by inspecting Eqs. (2) and (3). If 
q > p then C grows faster than B, so at high energies 
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For the choice q = 2{p — 1), the eigenvalue Ai is propor- 
tional to 1 / E and therefore plays the role of an effective 
mass term, even though no mass parameter is present at 
high energies. Note that imposing this choice requires 
the dominant coefficient in C to be CPT even. The null 
entries in the /j,-t block of h^ s and the fast-growing ee 
element guarantee maximal o v T mixing at high en- 
ergies, consistent with observations of atmospheric neu- 
trinos [57-59]. For a suitable choice of the ratio B 2 /C, 
as discussed in the next section, the seesaw mechanism 
also reproduces the L/E oscillation signature in the SK 
experiment [60]. 

Since the elements of h^ s are real, the probability 
P Vh ^ Va of oscillation from Vf, to v a can be written in the 
simple form 



p vb^va =$ab-4:^2 U a > a U a > b UvaU b > b sin 2 {A a , b ,L/ '2), 

a'>b' 



(7) 

where the quantities A a >b> = \ a ' — \' arc the eigen- 
value differences and L is the baseline. For each flavor 
pair a, b, the above sum contains three terms labeled 
by the values of a', b 1 < a'. Each term is the prod- 
uct of an amplitude —AUUUU with a sinusoidal phase. 
The antincutrino-oscillation probabilities Pp b -^v a are ob- 
tained by CPT conjugation. Since A, B, and C are real, 
all processes are T invariant. As a result, CP violation 
occurs if and only if CPT violation does. Notice that 
CP-violating effects can appear even though no analogue 
of the phase S in the 3^SM exists in the puma model. 

All the above properties are insensitive to the ee com- 
ponent of the B term in h v cS . As a result, a modified 
texture h'^ s can be constructed in which the ee entry 
in the B term vanishes. We have verified that most of 
the properties discussed in the remainder of this work 
remain unchanged for this modified texture. One excep- 
tion is the renormalizable model presented in Sec. VA, 
for which we use a zero ee entry in the B term because 
the nonzero value produces a tension between the de- 
scriptions of long-baseline reactor and of solar neutrinos. 



III. EXPERIMENTS 

Next, we study the implications of the general model 
(1) for different experiments. Many characteristics of the 
model are generic. For definiteness, in this section we 
illustrate the discussion with a specific c%a$m model [9]. 
Some comments on variant models are provided in Sec. 
IV. 

The numerical values of the three parameters in the 
Csa^m model are 

m 2 = 2.6 x 1(T 23 GeV 2 , 
a (5) = -2.5 x 1(T 19 GOV" 1 , 
1.0 x 10~ 16 GcV~ 4 , 



*(8) 



(8) 

The nonzero value of implies this model contains 
CPT violation. The value for m 2 is consistent with limits 
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FIG. 1: Energy dependences of the oscillation lengths for neu- 
trinos (top) and antineutrinos (bottom). The disappearance 
lengths for the puma model are L31 (top, solid line), L21 
(top, dashed line), L31 (bottom, solid line), and L21 (bottom, 
dashed line), displayed for the values (8). The dotted lines are 
the disappearance lengths L (solar) and L a tm (atmospheric) 
in the 3^SM. 




FIG. 2: Flavor content of the three neutrino eigenstates of 
ftg ff (left) and the three antineutrino eigenstates of h^ s (right) 
as a function of energy. For the puma model, the left-hand 
panel shows the energy dependences of | U a i e \ 2 (white), |t/ a ' M | 2 
(light grey), and |l7 a v| 2 (dark grey) for each neutrino mass 
eigenstate u a i , a 1 — 1,2,3, while the right-hand panel dis- 
plays the analogous energy dependences for antineutrinos. 
For the 3^SM, the corresponding quantities |C/ a / e | 2 (regions 
above dashed lines), |t/ a ' M | 2 (regions between dashed and solid 
lines), and |C/ a ' T | 2 (regions below solid lines) for neutrinos and 
those for antineutrinos are energy independent. The models 
coincide at all energies for the eigenstates vz, V3, but v\, V\ 
match U2, V2 only at low energies. 



General features 



from direct mass measurements and cosmological bounds 
[!]■ 

By construction, and c' 8 -* are the only nonzero 
SME coefficients denned in an isotropic frame /. In some 
scenarios, it is reasonable to identify I with a universal 
inertial frame U such as that denned by the cosmic mi- 
crowave background (CMB), but other possibilities exist. 
Whatever the choice for /, the experiment frame E is 
boosted in it by some combination of the Earth's motion 
relative to the CMB, the Earth's revolution about the 
Sun, and the Earth's rotation. The coefficients d' 5 ) and 
c^ 8 ) therefore induce anisotropic effects via the net boost 
in /. These could, for example, be detected by searches 
for sidereal or annual variations in E [42] . Experimental 
constraints and signals must be reported in a specified 
frame, but the frame E itself is inappropriate because 
it is noninertial and experiment specific. By convention, 
the canonical inertial frame used to report results is a 
Sun-centered frame S [15, 61]. Inspection reveals that 
the size of the effects in S induced by the values (8) all 
lie below the sensitivity levels achieved in experiments 
to date [36, 38-41]. Future experiments might offer im- 
proved sensitivity and thereby provide a distinct avenue 
for testing the model. 



The predictions of any model for neutrino and antineu- 
trino oscillations can be visualized using a certain plot in 
E-L space [8]. Experiments are represented on the plot 
as regions determined by their baseline and energy cover- 
age, while a given theory is represented by its characteris- 
tic oscillation wavelengths L a iy = 27r/|A a /(/| associated 
with the eigenvalue differences A. a iy{E). The absolute 
value is used because the oscillation phase is insensitive 
to the sign of A a 'b>. Each curve L a iy = L a n,i{E) indi- 
cates the first maximum of a kinematic phase in the oscil- 
lation probability, thereby establishing the minimal dis- 
tance from the neutrino source required for appearance 
or disappearance signals in a specific oscillation channel. 
Substantial signals appear in the region above each curve 
but are suppressed below it. 

Figure 1 shows this plot for the puma model with 
values (8) and the 3^SM. The 3^SM has two indepen- 
dent oscillation lengths, L Q — iirE/ArriQ and L atm = 
47T.E / Am^ tm , both of which grow linearly with the en- 
ergy and are therefore represented by straight lines in 
the plot. In the puma model, however, the unconven- 
tional energy dependences from B{E) and C(E) produce 
more general curves instead. These curves partially dif- 
fer for neutrinos and antineutrinos, a consequence of the 
CPT violation implied by the values (8). 
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The figure shows that the puma curves merge with 
the 3^SM lines Lq and L atm at low and high energies, 
respectively, suggesting consistency of the puma model 
with results in KamLAND, solar, and atmospheric ex- 
periments. This agreement is confirmed in the subsec- 
tions below. However, the two models are qualitatively 
different at intermediate energies. 

Novel effects arise from the unconventional energy de- 
pendence of /ig ff , which generates energy- dependent mix- 
ing. The flavor content of the three eigenstates of h^ s 
therefore changes with energy. Figure 2 shows this energy 
dependence for the values (8). At low energies, the flavor 
content approaches the tribimaximal limit. However, at 
high energies the eigenstate u 2 becomes completely popu- 
lated by v e . This implies the mixing <-> v T is maximal 
and controlled by A31 . The onset of this feature coincides 
with the onset of the Lorentz-violating seesaw. Indeed, 
as the mass term A becomes negligible in h" s , the frac- 
tion of v e in h>2 grows with the separation between the 
lines L21 and L31 in Fig. 1. 

Notice that the mixing angles in the 3^SM are en- 
ergy independent parameters that can freely be chosen 
to match data. In contrast, the mixing angles in the 
puma model at low and high energies are determined by 
the texture of h^ s and therefore are fixed features of the 
model that cannot be adjusted according to experiment. 
This reduced freedom is one reason why the puma model 
offers a more economical description of confirmed neu- 
trino data than the 3^SM. 

The energy dependence of the mixing matrix U implies 
the oscillation amplitudes —4UUUU in each flavor chan- 
nel and the corresponding probability (7) are also energy 
dependent. For given flavors a, b, the oscillation ampli- 
tudes are shown in Fig. 3. Note that negative amplitudes 
occur for disappearance channels, while positive ampli- 
tudes occur for appearance channels. The S 2 symmetry 
of /ig ff implies the four amplitudes for v e — > v T , v e — > V T , 
v T — > v T , and v T — > V T arc identical to those shown in the 
corresponding four central panels in the figure. The low- 
energy S3 symmetry of h" s forces the low-energy ampli- 
tudes to values set by tribimaximal mixing and ensures 
the low-energy equalities P v ^ Ve = = P Vt ^ Vt 

and Pv e ^v e = Pv^v^ = Pv T ^v T - At high energies, the 
amplitudes become either zero or one due to the Lorentz- 
violating seesaw mechanism. The lower four panels in the 
figure reveal that the dominant amplitude at high ener- 
gies has (a',b') = (3,1), leading to maximal <-» v T 
mixing and to an oscillation phase proportional to A31 
and hence to 1/E. Note also that the zero component 
U3 e of the mixing matrix (3), which is a consequence of 
the null eigenvalue of h^ s , implies that the oscillation in 
any channel involving v e or v e is controlled by only one 
amplitude because the other two vanish. 

The three figures reveal many of the evolving proper- 
ties associated with and h v eS at intermediate energies. 
For example, a peak appears between 10 MeV and 100 
MeV in the L 31 curve for neutrinos in Fig. 1, accompa- 
nied by corresponding features in Figs. 2 and 3. The 
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FIG. 3: Energy dependence of the oscillation amplitudes in 
the puma model. In each flavor channel, the amplitude factors 
-4UUUU in Eq. (7) are plotted for each of the three (a', 6') 
values 21 (solid lines), 31 (dashed lines), and 32 (dotted lines). 



peak represents a divergence in L31, which occurs when 
A a 'b' vanishes. Using the exact expressions (2) for the 
eigenvalues of h" s , we find that in general peaks occur 
for all positive energies E solving the equation 



A(B ~C) + B 2 



0. 



(9) 



The peaks can in general occur for both neutrinos and an- 
tineutrinos. The absence of these features in the antineu- 
trino plots suggests an origin in CPT violation. Since 
coefficients for CPT-odd Lorentz violation reverse sign 
under a CPT transformation, the nature of the solutions 
to Eq. (9) for antineutrinos changes. For the values (8), 
a single positive energy solves this equation for neutri- 
nos, but no solutions exist for antineutrinos and hence 
no antineutrino peaks arise in Fig. 1. 



B. Reactor antineutrinos 

In the puma model, the general survival probability for 
reactor antineutrinos is 

Pv^v e = 1 - 16 ( l_tf 2 )4 sin ' (5 A 2i^) • (10) 
N 1 N 2 

At low energies, the A term in h^ s dominates. Using 
the low-energy limits N ± — > (SA 2 , N 2 — > 3A 2 , we find for 
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FIG. 4: Reactor-antineutrino survival probabilities as a func- 
tion of L/E in the puma model (solid line) and in the 3^SM 
(dashed line). The data are from the long-baseline Kam- 
LAND experiment, for which L ~ 180 km [54]. 



Pv e ^v e the simple low-energy approximation 

Pv.-+v. « 1 - I sin 2 {^^j (low energy). (11) 

The fixed value 8/9 for the oscillation amplitude matches 
expectations because at low energies is diagonalizcd 
using the tribimaximal mixing matrix. This result ap- 
plies to reactor antineutrinos in both long- and short- 
baseline experiments. The large disappearance ampli- 
tude for reactor antineutrinos is evident in the 1> e — > V e 
panel of Fig. 3. 

1. Long-baseline reactor: KamLAND 

In the 3^SM, the reactor-antineutrino survival proba- 
bility for long-baseline experiments is 

J*££ « 1 - sin 2 29 12 sin 2 ■ (12) 

The data indicate values for the 3^SM parameters of 
sin 2 2<9 12 ~ 0.92 and Am| ~ 7.58 x 10" 5 eV 2 [53]. Com- 
paring the oscillation phase in this result with that in 
Eq. (11), we find that agreement with the KamLAND 
results can be achieved by choosing the mass parameter 
m 2 to be m 2 = Amg/3 [9]. This gives the numerical 
value adopted in Eq. (8). The match between the two 
models is shown in Fig. 4. 

Notice that the disappearance of reactor antineutrinos 
is described using only one parameter m instead of the 
usual two in the 3fSM. The conventional solar mixing 
angle 6 12 is eliminated as a degree of freedom by the 
form of the texture h^ s . Inspecting Eq. (11) reveals that 
at low energies the effective value of sin 2 29i 2 is numer- 
ically fixed to (sin 2 29i 2 ) c s ~ 0.89, which is close to the 
measured magnitude. The reader is however cautioned 
that this interpretation fails at higher energies due to the 
energy dependence of the mixing in the puma model. 



2. Short-baseline reactors 

In recent years, numerous reactor experiments with 
short baselines L < 1 km such as Bugey (L ~ 15,40 m) 
[62], CHOOZ (L ~~1 km) [63], Gosgen (L ~ 38, 46, 65 m) 
[64], and Palo Verde (L ~ 750,890 m) [65] have sought 
evidence for the disappearance of electron antineutrinos 
with null results. The explanations of these results differ 
qualitatively in the puma model and the 3^SM. 

In the 3z/SM, the z7 e survival probability is 

P^SM „ 1 sin 2 20i3 sin 2 (^|^ . (13) 

For energies E ~ 3 MeV, this gives an antineutrino disap- 
pearance length 27rS/Am 2 tm ~ 8.1 x 10 18 GeV -1 , which 
is about 1.5 km. The null experimental results are there- 
fore interpreted in the 3i^SM as a consequence of a small 
mixing angle 6*13. Note that the 3^SM survival prob- 
abilities (12) and (13) for long and short baselines, re- 
spectively, have the same form but involve four different 
parameters, Am 2 olar , 9 12 , Am 2 tm , and 9 13 . 

In contrast, in the puma model the oscillation proba- 
bility (11) holds at low energies for any baseline. Only 
the single parameter m is required to describe both the 
long- and short-baseline data. For energies E ~ 3 MeV, 
the antineutrino disappearance length is L 2 \ = 7r/A 2 i ~ 
27r£7Am| ~ 2.5 x 10 20 GeV" 1 , which is about 50 km. 
The null reactor results are therefore understood in this 
model as a consequence of the short baselines, which limit 
the contribution of the oscillation phase to the survival 
probability, rather than a consequence of a small oscilla- 
tion amplitude as in the 3^SM. Indeed, the amplitude of 
the oscillating term in Eq. (11) is 8/9, which is large. 

Since the puma model contains no term with a phase 
involving Am 2 tm , we see that at low energies the effec- 
tive value (sin 2 2#i3) e ff of the 3^SM quantity sin 2 2#i3 is 
exactly zero. This is a consequence of the zero value 
of Use, as can be confirmed by comparing the 3^SM 
mixing matrix with the tribimaximal limit of the mix- 
ing matrix (3). Note, however, that the energy depen- 
dence of the mixing matrix makes this result invalid at 
higher energies, where the effective value (sin 2 26* 13 ) c ff 
extracted from high-energy experiments can be nonzero 
even though Z7 3e identically vanishes. 



C. Solar neutrinos 

For neutrinos propagating in matter, the effective 
hamiltonian h" s acquires an additional term [66]. The 
modified effective hamiltonian {h" s ) M in the solar inte- 
rior can be written as 

(K S )ab - (KftU + V S ae S be , (14) 

where the solar matter potential Vq takes the value Vq = 
V2G F n e ~ 7.84 x 10~ 21 GeV at the solar core [67]. 
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FIG. 5: Averaged survival probability for solar neutrinos in 
the puma model (solid line) and in the 3^SM (dashed line). 
Both cases include matter-induced effects in the adiabatic ap- 
proximation. The data are from Ref. [52]. 
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FIG. 6: Survival probability for atmospheric neutrinos as a 
function of L/E in the puma model (solid line) and in the 
3^SM (dashed line) compared to SK data, for which Lsk — 
600 km [60]. 



The presence of the solar potential preserves the puma 
texture (1) because it corresponds to a simple redefinition 
of the function C of the form C — >■ C + Vq . The exact 
eigenvalues and the exact mixing matrix in the presence 
of matter can therefore be found immediately by apply- 
ing this redefinition to Eqs. (2), (3), and (4). A short 
calculation reveals that the averaged survival probability 
of solar neutrinos takes the exact form 



(Pv e -H> e ) 



(Af - 2A) (Ai - 2A) 



N M 



(Xf - 2A) (A 2 - 2A) 



No 



(15) 



For the lower-energy region of the solar spectrum with 
E <~ 0.1 MeV, the solar potential V Q and the functions B 
and C are negligible. In this limit, the averaged survival 
probability becomes 



(Pv e ^v e ) « ^2\U a > e \ 4 = § (low energy) 



(16) 



in agreement with the data. This result is to be expected 
because the vacuum mixing matrix is tribimaximal at low 
energies. 

For higher energies, the solar potential and the 
Lorentz-violating terms can introduce novel effects, de- 
pending on the form of the functions B and C. The de- 
tailed form of the averaged survival probability therefore 
becomes model dependent. However, the neutrino sur- 
vival probability initially drops below the limiting value 
5/9 as the energy increases. This generic effect is a con- 
sequence of the energy independence of V©, which en- 
sures Vq becomes relevant at energies comparable or be- 
low those for the Lorentz-violating terms and thereby 
enhances the disappearance of v e . 

The above features are visible in Fig. 5. The solid line 
displays the averaged survival probability for the values 
(8). The curve is similar to that obtained from the 3fSM 
and is compatible with observations. 



D. Atmospheric neutrinos 

In the puma model, the exact survival probability of 
atmospheric neutrinos is 



1 -4 

-2 
-2 



(A + B) A 
N 2 N 2 
(A + B) 2 . ., , 

(A + B) 
N 2 



sin 2 (±A 21 L) 



in 2 (±A 31 L) 



sin 



2/1 



A 32 £) • 



(17) 



However, as E grows the A term becomes negligible, so 
the sole mass parameter m is irrelevant for high-energy 
oscillations. Requiring C to increase with energy faster 
than B yields the high-energy limits Ai — > —2B 2 /C, 
A 2 -> C, N 2 -> 2B 2 , and N% -> C 2 . The electron- 
neutrino content then lies exclusively in the second eigen- 
state. This leaves the other two uniformly populated by 
and v T , as can be verified by examining Fig. 2. The 
survival probability (17) then takes the simple form 



1 - sin" 



B 2 L 



C 



(high energy) 



(18) 



where in the second equation the leading contributions 
to B and C are expressed as monomials of the form (5). 
Note that the unit amplitude of the oscillation term im- 
plies maximal mixing, as discussed following Eq. (6). 

In contrast, the 3^SM survival probability for atmo- 
spheric neutrinos takes the form 



P. 



~ 1 - sin 2 26» 23 sin 2 



4E 



(19) 



depending on two parameters # 2 3 
imental data provide the values sin 2 2# 23 > 0.90 and 



and Am 2 tm . Exper- 



8 



\Aml tm \ ~ 2.32 x 10~ 3 cV 2 [57]. Comparison of Eqs. 
(18) and (19) suggests agreement with atmospheric data 
can be obtained when the ratio of (k^) 2 and satisfies 

This condition has been used to constrain the coefficients 
and in Eq. (8) [9]. The resulting match between 
the two models is shown in Fig. 6 along with SK data. 

Note that the ratio (20) represents only one degree of 
freedom. Nonetheless, it suffices to reproduce the data 
for atmospheric neutrinos via Eq. (18). The other degree 
of freedom in the two coefficients k^ p \ determines 
the onset of the Lorentz- violating seesaw. Increasing 
while holding fixed the ratio (20) causes the seesaw to 
trigger at lower energies. 
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FIG. 7: Comparison of the puma model (solid lines; xl = 1-0, 
Xtt = 0.9), the tandem model [19] (dotted lines; xt = 1-9, 
xl = 1.0), and the 3i/SM (dashed lines; xl = 2-2, xl = 1-1), 
with MiniBooNE neutrino [5] and antineutrino [6] data. 



E. Short-baseline accelerator neutrinos 

At high energies E > 1 GeV, a variety of short-baseline 
experiments have reported null results. BNL-E776 (L = 
1 km) searched for — > v e and 77^ — > V e at 1 GeV 
[68]. CCFR (L ~ 1 km) searched for -> v e , V ^ -> 
77 e , v e -> z/,-, and 77 e -> z7 T at 140 GeV [69]. CDHS 
[L ~ 130 m) searched for disappearance at 1 GeV 
[70]. CHORUS (L ~ 600 m) searched for ^ -> z/ r at 27 
GeV [71]. NOMAD (L ~ 600 m) searched for z^ — > i/ T 
and i/ e -> z/ T at 45 GeV [72]. NuTeV (L ~ 1 km) searched 
for z/^ -4- f e and 77^ -> z7 e at 150 GeV [73]. 

The puma model is consistent with all these null re- 
sults. For energies above the seesaw scale ~ 1 GeV, 
-H> v T mixing becomes maximal by construction, as 
described following Eq. (6). This feature implies vanish- 
ing high-energy mixing and hence no oscillations in the 
channels — > v F , — > V e , v e — > v T , and v e — > v r . The 
behavior can be seen directly from Fig. 3, which displays 
the energy dependence of the oscillation amplitudes. 

In the Vfj, — > v T channel, the oscillation amplitude 
is maximal at high energies. However, the oscillation 
phase is controlled by A 2 i, which generates an appear- 
ance length L21 of several hundred kilometers at 1 GeV. 
The lack of a signal in this channel in the CHORUS or 
NOMAD data is therefore understood here as a conse- 
quence of their short baselines. 



F. MiniBooNE anomalies 

Two results from the MiniBooNE experiment indicate 
possible oscillation effects that cannot be accommodated 
within the 3^SM. For neutrino oscillations -4- v e , 
MiniBooNE finds a 3<t excess of events at low energies 
around 200-500 MeV [5]. For antineutrino oscillations 
F p — > V e , a 1.3cr low-energy excess has also been reported 
[6], with recent preliminary data suggesting a larger ex- 
cess [74]. 



These results are interesting in the present context 
because they lie in the energy region where the seesaw 
mechanism is triggered. Following onset of the seesaw, 
the eigenvalue Ai decreases linearly with energy while A2 
grows rapidly. The appearance length L21 oc (A2 — Ai) -1 
therefore drops steeply, becoming a few hundred meters 
at MiniBooNE energies. This produces a large oscillation 
phase and hence a signal in the experiment. However, the 
oscillation amplitude for o v e mixing rapidly goes to 
zero as the -h- v t mixing becomes maximal, as can 
be seen in Fig. 3. As a result, the appearance signal in 
MiniBooNE vanishes at higher energies. 

The puma model therefore naturally describes a low- 
energy excess in MiniBooNE. Moreover, the excess can 
differ substantially for neutrinos and antineutrinos when 
a coefficient for CPT-odd Lorentz violation is involved, 
such as occurs for the c%a^m example (8). In general, 
the energy at which the excess appears depends on the 
seesaw scale and becomes smaller as increases. For 
the values (8), the match to data is shown in Fig. 7. 

We emphasize that these interesting features of the 
model arise without introducing additional particles or 
forces. They are a consequence of the comparatively el- 
egant texture (1) that describes all compelling neutrino- 
oscillation data. 



IV. PREDICTIONS 

The discussion in the previous sections demonstrates 
that two of the three parameters of the model (1) suf- 
fice to reproduce all the compelling data for neutrino 
and antineutrino oscillations, while the third accommo- 
dates the two MiniBooNE anomalies. Comparison to the 
3^SM, which uses five nonzero parameters to describe es- 
tablished results but cannot reproduce the MiniBooNE 
anomalies, suggests the puma model offers a frugal inter- 
pretation of known data. 

The model predicts a variety of signals, some of which 
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FIG. 8: Prediction for the probabilities of v e appearance (left) 
and of ~D e appearance (right) in various long-baseline exper- 
iments according to the puma model (solid lines) and the 
3^SM (upper dashed lines, sin 2 2#i3 = 0.02; middle dashed 
lines, sin 2 2^13 = 0.05; lower dashed lines, sin 2 2#i3 = 0.08). 
Matter effects are included. 




FIG. 9: Prediction for the probabilities of f M disappearance 
(left) and of F M disappearance (right) in various long-baseline 
experiments according to the puma model (solid lines) and 
the 3z^SM (dashed lines). Matter effects are included. The 
data are taken from Refs. [7, 57, 79]. 



differ qualitatively from 3f SM expectations. In this sec- 
tion, we address some features of relevance to future ex- 
periments. 

A. Long-baseline neutrinos 

A variety of long-baseline experiments, including 
LBNE (L ~ 1300 km) [75], MINOS (L ~ 735 km) [76], 
NOi/A (L ~ 810 km) [77], and T2K (L ~ 298 km) [78], 
have design capabilites to search for v e appearance in a 
beam. These searches are motivated in part by the 
prospects of measuring the 3^SM parameter 6*13. No such 
parameter exists in the model (1), but signals in these ex- 
periments may nonetheless appear. 

To characterize potential signals in the puma model, 
recall that the appearance length L21 for — > v e de- 
creases steeply with energy due to the seesaw mechanism. 
Baselines L >• L 2 i therefore involve rapid oscillations, so 
accelerator experiments with long baselines can observe 
only the averaged oscillation probability, given exactly 

by 

(P _ )-l {A + BY (21) 

To allow for matter effects on neutrinos traversing the 
Earth, Nf 4 and are given by Eq. (4) with the re- 



placement C — > C + V®, where the Earth's matter po- 
tential Vfs is Vfs ~ 1.2 x 10~ 22 GeV. 

For energies above the seesaw scale, -B- v T mixing 
dominates while -R- v e mixing is highly suppressed. 
However, to describe the SK and MINOS data, the see- 
saw must trigger below 1 GeV. This means only small 
signals from — > v e transitions can appear in the high- 
energy experiments LBNE, MINOS, and NOi^A. In con- 
trast, T2K runs at lower energies, and so a larger appear- 
ance signal that decreases rapidly with the energy is to 
be expected. Quantitative predictions for the probabili- 
ties for v e and z7 e apppearance in the various experiments 
are shown in Fig. 8 for the values (8). Note that matter 
effects are almost negligible compared to the large eigen- 
value A2 controlling the mixing, whereas for the 3^SM 
curves they induce substantial differences between the 
probabilities for v e and v e appearance. 

We remark in passing that attempting to interpret 
these signals as arising from a nonzero 3z^SM angle #13 
would predict that #13 is larger in T2K than in the other 
higher-energy experiments. This is compatible with re- 
cent results for v e appearance [2]. Note also that within 
this perspective the effective values of 6*13 obtained with 
long-baseline accelerators are unrelated to the effective 
values of 6*13 extracted from studies of reactor antincu- 
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FIG. 10: Predictions for the probability of 77 e disappearance 
in the Daya Bay, Double Chooz, and RENO experiments ac- 
cording to the puma model (solid lines) and the 3^SM (dashed 
lines, labeled with the value of sin 2 2^13). 
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FIG. 11: Predictions for the probability of 77 e appearance 
at the three baselines proposed for the DAE5ALUS experi- 
ment according to the puma model (solid lines) and the 3^SM 
(dashed lines, labeled with the value of sin 2 2^13). 



trinos discussed in Sec. IIIB 2. 

High-energy long-baseline experiments can also per- 
form precision studies of disappearance. For this os- 
cillation channel, no differences between the puma model 
and the 3^SM are expected in LBNE, MINOS, and 
NOf A because they operate at energies above the seesaw 
scale. However, significant differences between neutrinos 
and antincutrinos are predicted for the lower-energy por- 
tion of the T2K spectrum when the theory contains a co- 
efficient for CPT-odd Lorentz violation, as in the Cga^m 
model. The predictions are displayed in Fig. 9, where as 
before the values (8) are used for illustration. 



B. Short-baseline neutrinos 

Using baselines in the range 1-2 km, modern reactor 
experiments such as Daya Bay (L ~ 1985 m) [80], Dou- 
ble Chooz (L ~ 1050 m) [81], and RENO (L ~ 1380 
m) [82] propose to measure the disappearance of reactor 
antineutrinos. Like their long-baseline cousins, these ex- 
periments are driven partly by prospects for measuring 
the 3^SM mixing angle 6*13. 

In the puma model, the T> e survival probability for en- 
ergies 2-9 MeV has a large amplitude, as can be seen 
from Eq. (11). The oscillation signal in a given reac- 
tor experiment therefore depends only on the size of the 
baseline L compared to the antineutrino disappearance 
length L21 — 50 km discussed in Sec. IIIB 2. The base- 
lines for the Daya Bay, Double Chooz, and RENO ex- 
periments are all short compared to this, with Daya Bay 
having the greater sensitivity to oscillation signals due 
to its longer baseline. Since L 2 i grows linearly with the 
energy, any oscillation signal in these experiments is ex- 
pected to appear predominantly at low energies. Using 



the model values (8), the predictions for the disappear- 
ance probabilities in the three experiments are shown in 
Fig. 10. 

We note in passing that the recent suggestion of an 
ovcrestimation of antineutrino fluxes in short-baseline re- 
actor experiments [83] is difficult to reconcile with the 
three-parameter puma model. Since effects at low ener- 
gies are governed by only one parameter m, which is fixed 
by KamLAND data, no other oscillation length appears 
at reactor energies. The existence of only one parameter 
is a consequence of the S3 flavor symmetry, so a slight 
breaking of this symmetry at low energies could accom- 
modate an additional parameter and hence a correspond- 
ing signal. This construction would introduce another 
degree of freedom but requires no additional neutrinos. 
However, investigations along these lines lie beyond the 
scope of the present work. 

Another experiment of interest in the context of short- 
baseline neutrinos is the recent DAE<5ALUS proposal [84] 
to study CP violation, which would generate neutrinos 
at several different baseline distances from a detector us- 
ing high-power accelerator modules to beam protons onto 
graphite sources. A popular configuration would offer the 
capability to search for F M — > v e transitions using three 
baselines of about 1.5 km, 8 km, and 20 km. The large 
oscillation amplitude in this region suggests appearance 
signals in the detector can be expected from the more 
distant sources. The predicted appearance probabilities 
for the three proposed baselines are shown in Fig. 11. 
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FIG. 12: Averaged survival probability for solar neutrinos in 
the CAd-ifn model (solid line) and in the 3^SM (dashed line). 
Both cases include matter-induced effects in the adiabatic ap- 
proximation. The data are taken from Ref. [52]. 



V. VARIANT PUMA MODELS 

In the preceding sections, the implications of the gen- 
eral texture (1) have been illustrated with the c 8 a 5 m 
model, using the specific values (8). However, other 
models can be constructed using h^ s that successfully 
describe most or all compelling neutrino data. Some of 
these offer distinctive features or intriguing possibilities 
for describing experimental anomalies beyond the Mini- 
BooNE ones. This section outlines some results for a few 
of these variant models. 




0.1 









^^i^v^/, ------ _ 




0.1 



0.2 0.4 0.6 

E (GeV) 



0.2 



0.4 0.6 
E (GeV) 



0.8 



FIG. 13: Prediction for the probabilities of v e appearance 
(left) and of V e appearance (right) in various long-baseline 
experiments according to variant puma models (solid lines, 
labeled by the model) and the 3^SM (upper dashed lines, 
sin 2 26<i3 = 0.15; lower dashed lines, sin 2 20 i3 = 0.05). Matter 
effects are included. 



A. The C4a-$m model 

The texture h^ s with the smallest monomial orders p 
and q in Eq. (5) requires only renormalizable operators 
of dimensions 3 and 4 in the minimal SME, hence pro- 
ducing a C4a 3 m model. For definiteness, we adopt in this 
subsection the specific numerical values 

m 2 = 2.6 x 1CT 23 GeV 2 , 
a (3) = -2.5 x 1(T 21 GeV, 
c (4) = 1.0 x 1(T 20 . (22) 

As mentioned at the end of Sec. II, we use a zero ee en- 
try in the B term for this model. The coefficient 
comes with CPT violation, so differences between neu- 
trino and antineutrino properties can be expected. These 
values are consistent with limits from direct mass mea- 
surements, cosmological mass bounds, and constraints on 
anisotropic oscillations. 

The C4a 3 m model is compatible with all accepted ex- 
perimental oscillation results discussed in Sec. Ill, includ- 
ing those obtained with reactor, solar, and atmospheric 
neutrinos. However, the eigenvalue A 2 grows too slowly 
to produce a signal in MiniBooNE because the function 
C is linear in energy. 



An interesting qualitative difference introduced by the 
model appears in the predicted averaged survival proba- 
bility for solar neutrinos, shown in Fig. 12. The probabil- 
ity curve incorporates a striking neutrino-disappearance 
maximum in the central-energy region, despite passing 
though all data points. This reflects the importance at 
lower energies of the coefficient c^ 3 ' , an effect absent for 
the coefficient in the csa^m model. The curve shape 
suggests future analyses of solar data in the 1-10 MeV 
part of the neutrino spectrum could provide an interest- 
ing experimental test of the model. 

Another distinctive feature of the model is a large sig- 
nal for v ^ — > v e oscillations in long-baseline experiments. 
The signal decreases with energy, as shown in Fig. 13. 
Analysis of the recent data supporting electron-neutrino 
appearance in the T2K and MINOS experiments [2] could 
provide a sharp constraint on this signal, potentially ex- 
cluding the values (22). 



B. The C6C4m model 

Other interesting variant models can be constructed 
using only CPT-even operators. Using monomials of 
smallest order produces a c^c^ra model. We choose here 
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FIG. 14: Comparison of the cqCatti model (solid lines; xi ~ 
1-0, X77 = 0.8), the tandem model [19] (dotted lines; xi = 1-9, 
X% = 1.0), and the 3^SM (dashed lines; xl = 2.2, x% = 1-1), 
with MiniBooNE neutrino [5] and antineutrino [6] data. 



the specific numerical values 



to 2 = 2.6 x 1CT 23 GeV 2 

c< 4 > = 7.7 x lCT 20 , 

c (6) = 1.0 x 1(T 17 GeV~ 



(23) 



As for other examples considered in this work, these 
values are consistent with limits from direct mass mea- 
surements, cosmological mass bounds, and constraints on 
anisotropic oscillations. 

Like the c%a^m and c^a^m models, the cq 04.712 model 
provides a good match to the data from reactor, solar, 
and atmospheric neutrinos discussed in Sec. III. The 
presence of the c^' term substantially affects the physics 
in the region 10 MeV to 1 GeV. It generates a signal in 
the MiniBooNE region that includes low-energy excesses 
in both neutrinos and antineutrinos, as shown in Fig. 14. 
Note the asymmetry in the signal, which here reflects 
experimental acceptance rather than CPT violation. 

For v e appearance in long-baseline experiments, the 
probabilities are generically closer in magnitude to those 
of the 3^SM with a moderate value of 6*13, as shown 
in Fig. 13. One interesting feature is the substantially 
larger signal produced in T2K relative to MINOS, which 
is compatible with the central values of recently reported 
measurements [2]. 



Four-coefficient models 



Some indication of a potential difference between 
and disappearance probabilities has recently been re- 
ported by the MINOS collaboration [7]. Although the ef- 
fect may disappear with improved statistics, a difference 
of this kind is of interest in the present context because 
it cannot be accommodated in the 3fSM, which requires 
identical neutrino and antineutrino masses. 

In three-parameter puma models, the survival proba- 
bilities of muon neutrinos and antineutrinos at MINOS 
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FIG. 15: Prediction for the probabilities of disappearance 
(left) and of F M disappearance (right) in various long-baseline 
experiments according to the c%a$cmi model (solid lines; xt = 
1.4, xl = 0.9 for MINOS) and the 3^SM (dashed lines; X l = 
1.0, xl = 1-6 for MINOS). Matter effects are included. The 
data arc taken from Refs. [7, 57, 79]. 



energies are generically the same. The result holds be- 
cause at high energies the relevant eigenvalue difference 
is A31 w 2B 2 /C. This is even under CPT provided the 
B term is either odd or even, which is true whenever the 
B term is only a monomial in E. For the cga^m model 
for example, this symmetry is reflected in Fig. 9. These 
puma models therefore cannot accommodate anomalies 
of the MINOS type either. 

In this subsection, we show that an ad hoc modifica- 
tion using an additional parameter can describe anoma- 
lies of this type. The key idea is as follows. Instead 
of choosing a monomial in energy for the B term, we 
can take a binomial involving two different monomial or- 
ders p and r < p, one even and one odd. This produces 
a four-coefficient model with both CPT-odd and CPT- 
even terms contributing at high energies. If the value of 
r is close to p and the corresponding coefficients are simi- 
lar in size, then the oscillation probabilities for neutrinos 
and antineutrinos differ at high energies. In particular, 
the energies of the first oscillation maxima of neutrinos 
and antineutrinos differ, which is a feature of the MINOS 
effect. However, to preserve compatibility with other ex- 
periments, the ee entry of the effective hamiltonian 
must remain unchanged. The extra coefficient should 
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therefore appear only in the e^i and er entries of h" s . 
One way to achieve this is to choose a binomial for the 
C term as well, compensating for the modification of the 
ee entry arising from the change to B. 

As an example, we can add a fourth coefficient to the 
Cga^m model while leaving unchanged its main features. 
Choosing r = 4, which satisfies the requirement r < p = 
5 with r near p, the fourth coefficient can be denoted 
c*- 4 -*. To generate a neutrino-antineutrino difference at 
high energies, it suffices to redefine the B and C terms 
as B -> B + c^E, C -> C - c^E. For definiteness, we 
can take the numerical values 



m — 
a< 6 > = 

*(8) = 



2.6 x 1(T 23 GeV 2 , 
2.0 x 1(T 20 , 
-2.6 x 1CT 19 GeV" 1 , 
1.0 x 10~ 16 GcV~ 4 , 



(24) 



which as before are consistent with limits from direct 
mass measurements, cosmological mass bounds, and con- 
straints on anisotropic oscillations. These choices pre- 
serve the attractive features of the simpler c^a^m model. 

The introduction of the fourth coefficient primarily af- 
fects the probabilities of and disappearance in long- 
baseline experiments, as shown in Fig. 15 for the values 
(24). In particular, the MINOS data in Rcf. [7] are rea- 
sonably described by the model. The location of the pri- 
mary minimum for antincutrino oscillations is at a higher 
energy than that for neutrinos, in agreement with the re- 
ported effect. We emphasize that this result is achieved 
with a single additional parameter, without any masses, 
in a global model of neutrino oscillations. Note also that 
this c^a^c^m model predicts a large difference between 
the probabilities of and disappearance in the T2K 
experiment. However, the v e appearance probabilities in 
long-baseline experiments are essentially unchanged from 
those for the c^a^m model shown in Fig. 8. 
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FIG. 16: Energy dependences of the oscillation lengths for 
antineutrinos in the doubly enhanced puma model. The cor- 
responding plot for neutrinos in the top panel of Fig. 1 re- 
mains unaffected. The disappearance lengths for the model 
are L31 (solid line) and L21 (dashed line), displayed for the 
values given by Eqs. (24), (27), and (28). The dotted lines are 
the disappearance lengths L (solar) and L a tm (atmospheric) 
in the 3^SM. 
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FIG. 17: Comparison of enhanced puma model (solid line, 
X 2 = 1.6) and the 3//SM (dashed line, X 2 = 2.6) with LSND 
antineutrino data [5]. 



D. Enhanced models 

In a search for appearance of electron antineutrinos 
in a beam of muon antineutrinos, the LSND experiment 
found evidence for a small probability P v ^ Ve ~ 0.26 ± 
0.08% of 77 M — > V e oscillations at baseline L = 30 m and 
energies in the range 20-60 MeV [4]. This signal cannot 
be accommodated within the 3^SM because the required 
mass-squared difference Atol SND ~ 1 eV 2 is orders of 
magnitude larger than Amg and Ara 2 tm . 

In the general puma model, the oscillation v ^ — > 77 e is 
given exactly by 

(25) 

N 1 N 2 \ £ / 

which is governed by the same eigenvalues that control 
reactor antineutrinos. The energies in the LSND experi- 
ment are greater than those of reactor antineutrinos, but 



the appearance length is smaller by about two orders of 
magnitude. Achieving this with a monomial energy de- 
pendence in /i£ ff while preserving consistency with other 
experiments is challenging, as it requires a large power of 
the energy and a seesaw triggered around 10 MeV. 

An interesting option generating the required steep fall 
and rise in Li\{E) is to introduce a smooth nonpolyno- 
mial function with a peak. A function of this type could 
act as an enhancement of h" s arising from a series of co- 
efficients in the SME. It can be approximated generically 
using at least three parameters, one to position it, one 
to fix its height, and one to specify its width. A simple 
example is a gaussian enhancement of the form 

5h = aexp[-l3(E-e) 2 ]. (26) 

To preserve the £2 symmetry of h^ s , the enhancement 
can be limited to the e~p and er entries of h^ s via the 
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FIG. 18: Predictions for the probability of V e appearance at 
the three baselines proposed for the DAEJALUS experiment 
according to the enhanced puma model (solid lines) and the 
3f SM (upper dashed lines, sin 2 2#i3 = 0.15; lower dashed 
lines, sin 2 26»i3 = 0.05). 
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FIG. 19: Prediction for the probabilities of disappearance 
(left) and of disappearance (right) in various long-baseline 
experiments according to the doubly enhanced puma model 
(solid lines; X l = 1-3, xl = 0.9 for MINOS) and the 3vSM 
(dashed lines; xl = 1-4, X% = 1-8 for MINOS). Matter effects 
are included. The data are taken from Refs. [7, 57, 79]. 



redefinitions B — > B + Sh and C — > C — Sh. Under the 
CPT transformation mapping to h^ s , the signs of 
a and e change. As a result, either the function or its 
CPT conjugate is localized at an unphysical value of the 
energy, and so the enhancement affects either neutrinos 
or antineutrinos but not both. 

As an example, consider an enhanced cga^m model 
with specific enhancement values 

a = 3.0 x 10~ 19 GeV, 
P = 3.0 x 10 3 GcV~ 2 , 
e = 60 MeV. (27) 

The positive value of e ensures that only antineutrinos are 
affected. This enhancement produces a sharp dip in the 
disappearance length L,2\(E) centered around 60 MeV, as 
displayed in Fig. 16. The resulting oscillation probability 
Pv^v c contains a nonzero signal in the same region as 
the LSND data, as shown in Fig. 17. Since the L21 curve 
passes through the region of sensivity for the DAE<5ALUS 
experiment [84], large signals are predicted in all three 
detectors as shown in Fig. 18. A large oscillation signal is 
also predicted at high energies in the OscSNS experiment 
[85]. 

Differences between and disappearance proba- 
bilities can also be generated by another enhancement of 
this general type. Two simultaneous enhancements can 
be included without interference provided they are local- 
ized in different regions of the spectrum. Figure 16 shows 
the effect on the disappearance length L31 of adding a 
second enhancement with the values 

a 2 = -2.0 x 10~ 19 GeV, 
h = 13 GeV" 2 , 

e 2 = 1.7 GeV. (28) 

The resulting disappearance probabilities for muon neu- 
trinos and antineutrinos in long-baseline experiments are 
shown in Fig. 19. While the single extra coefficient in 
the four-parameter c^a^c^m model is more economical 
in generating an anomaly like that reported by MINOS 
[7] , the introduction of the second enhancement centered 
near 2 GeV produces interesting and distinctive oscilla- 
tion signals in the LBNE and NO^A experiments for 
disappearance, as shown in Fig. 19. In contrast, the v & 
appearance probabilities displayed for the c%a§m model 
in Fig. 8 are largely unaffected by the enhancement. 



VI. DISCUSSION 

In this work, we have investigated the behavior of neu- 
trinos governed by an effective hamiltonian h^ s of the 
puma form (1). This texture is interesting in part be- 
cause it leads to descriptions of neutrino oscillations that 
are globally compatible with experimental data. The as- 
sociated Lorentz- violating models are intriguing because 
they are frugal, they have a certain elegance, and their 
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novel features are compatible with data in unexpected 
ways. We remark in passing that the existence of these 
models was unclear a priori, becoming apparent only 
through a systematic search for viable candidates. 

The frugality can be traced to the use of only two de- 
grees of freedom to describe established data, instead of 
the usual five required by the 3^SM, while a third de- 
gree of freedom efficiently encompasses the MiniBooNE 
anomalies. Adding a fourth degree of freedom readily 
generates an anomaly of the MINOS type, while a three- 
parameter enhancement produces a signal in the LSND 
experiment. These four latter degrees of freedom are ad 
hoc, and their necessity depends on the ultimate confir- 
mation of the MINOS and LSND anomalies. However, 
to our knowledge the resulting texture represents the sole 
extant global model for neutrino oscillations, and more- 
over uses degrees of freedom comparable in number to 
those of the 3vSM. 

The symmetry of h^ s also implies a certain elegance. 
The puma texture (1) could naturally stem from more 
fundamental physics at the unification scale that gen- 
erates a democratic contribution to the dominant mass 
operator in the low-energy effective theory. The result- 
ing S3 symmetry then holds at low energies and ensures 
tribimaximal mixing but is broken to S2 at higher ener- 
gies by subdominant terms in the SME. This symmetry 
structure leads to the attractive quadratic calculability 
of the models. The coefficients required for compatibil- 
ity with data are of plausible Planck-suppressed size. 

The novel features of the puma models originate in the 
unconventional energy dependence in the eigenvalues of 
/ig ff and the mixing matrix U. Indeed, it is a pleasant sur- 
prise that the models pass the test of compatibility with 
existing data, despite their qualitative differences com- 
pared to the 3;/SM. One striking feature is the Lorentz- 
violating seesaw, which makes viable the absence of a 
mass parameter at high energies. Another satisfying fea- 
ture is the steep drop with energy of the oscillation length 
L21, which is naturally enforced by the third degree of 
freedom required to generate the Lorentz-violating see- 



saw. As discussed above, this drop enables L21 to at- 
tain the vicinity of the MiniBooNE experiment in E-L 
space, thereby generating a low-energy signal compatible 
with the MiniBooNE anomaly. Moreover, this feature 
appears in conjunction with a rapid decrease in the rele- 
vant oscillation amplitudes accompanying the large oscil- 
lation phase. This accounts for null signals in high-energy 
short-baseline accelerator experiments without invoking 
the tiny oscillation phase of the 3fSM. 

The puma texture h" s predicts certain signals that dif- 
ferentiate sharply between it and the 3^SM. One key fea- 
ture is the energy dependence of the effective mixing an- 
gle #13. This implies the probability of v e appearance is 
larger in the T2K experiment than in the MINOS exper- 
iment. It also predicts no accompanying signal in reactor 
experiments, a result at odds with the 3fSM. With an en- 
hancement present, strong signals are predicted in exper- 
iments at intermediate energies and moderate baselines 
such as the proposed DAE<5ALUS experiment. Another 
unique signal predicted by some models is CPT viola- 
tion, which implies differences in oscillation probabilities 
between neutrinos and antineutrinos. Perhaps the most 
direct evidence for Lorentz violation would be the dis- 
covery of oscillation anisotropics arising from the boost 
relative to the isotropic frame. One signal would be side- 
real variations of oscillations in the laboratory frame [42] , 
which in the puma models are predicted to be some 10- 
100 times below current sensitivities [9]. In any event, 
the results in this work show that Lorentz- and CPT- 
violating models can serve as an experimentally viable 
foil to the 3z/SM, while offering a simple and credible 
alternative for realistic modeling of neutrino oscillations. 
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